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ABSTRACT 

We present the high-resolution Chandra X-ray Observatory persistent (non-dip) spectrum of 
4U 1916—05 which revealed narrow absorption lines from hydrogenic neon, magnesium, silicon, and 
sulfur, in addition to the previous identified hydrogenic and helium-like iron absorption lines. This 
makes 4U 1916—05 only the second of the classical X-ray dipper systems to show narrow absorption 
lines from elements other than iron. We propose two possible explanations for the small measured 
line widths (< 500-2000 km s""'^), compared to the expected Keplerian velocities (> 1000 km s^^) of 
the accretion disk in this 50-min orbital period system, and lack of wavelength shifts (< 250 km s~^). 
First, the ionized absorber may be stationary. Alternatively, the line properties may measure the rel- 
ative size of the emission region. From this hypothesis, we find that the emission region is constrained 
to be ^0.25 times the radial extent of the absorber. Our results also imply that the ionized absorber 
spans a range of ionization parameters. 

Subject headings: binaries: close — stars: individual (4U 1916—05) — X-rays: binaries 
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1. INTRODUCTION 

It was anticipated that high-resolution spectra of low- 
mass X-ray binaries (LMXBs) would reveal a wealth of 
line features which could be used to better understand 
accretion disk structure, but such features have been 
largely absent in most Chandra X-ray Observatory and 
XMM-Newton spectra of LMXBs. One exception is the 
class of X-ray dippers. X-ray dippers show intensity dips, 
although not necessarily eclipses, on the orbital period 
and have inclination angles in the range 60°<i<80°. The 
dipping behavior is caused by absorption of the contin- 
uum by the accretion disk or its atmosphere. 

The best studied of the X-ray dippers, the neutron star 
(NS) binary, EXO 0748-676 (Porb=3.8 hr), showed emis- 
sion lines from the hydrogenic and helium-like species 
of nitrogen, oxygen, neon, magne sium, and sili c on in 
its XMM and Ch andra spectra l|Cottam et all 1200 It 
iJimenez-Garate et al. 2003). The XMM spectrum of the 
X-ray dipper MXB 1659-298 (Porb=7.1 hr) showed nar- 
row absorption lines from highly ionized o xygen, neon, 
and i ron, and a broad Ye-K emission line (|Sidoli et al.l 
120011) . XMM data from four other dippers, X1624-490 
(Porb=21.0 hr), 4U 1254-69 (Porb=3.9 hr), 4U 1916-05 
(P„,b=0.83 hr), and 4U 1323-62 (Porb=2.94 hr) re- 
vealed similar narrow iron absorption li nes and some 
evidence for broad iron emi s sion lines iParmar et al.l 
mV^ IRnirin fc PaTmaHl2003l IRoirin et al .112004 I200,5t 
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The presence of these absorption features in X-ray dip- 
pers, but not in other LMXBs, has led to the suggestion 
that the absorbing material has a cylindrical distribu- 
tion fe.g.. iBoirin et al1l2004|) . The lack of variation in 
the properties of the absorption lines as a function of 
orbital phase, excluding the dipping region, implies an 
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azimuthal symmetry (see, e.g.. IBoirin et all 12004'. and 
references therein). This in turn points to some rela- 
tionship between the accretion disk and the absorbing 
material. 

In this paper, we present the high-resolution Chandra 
spectrum of 4U 1916—05, the shortest orbital period X- 
ray dipper. The 50-min orbit al period of 4U 1916—05 was 
discovered independ ently bv lWhite fc Swan k (1982) and 
iWalter et alJ l)1982() from dips seen in the X-ray emission. 
The stability of the dip period over many years led the 
authors to associat e it with t he orbital period of the sys- 
tem (but also see , iChou et al. 2001: Homer et al. 2001t 
iRetter et alEoO^ . The dipping behavior of 4U 1916-05 
is variable on the order of ~4 days, varying from short 
(wlO% of orbital phase), shallow d ips to long (r:!40% 
of orbital phase), deep dips (e.g., iChou et aLl l200l(l . 
Unique among the X-ray dippers, the 50-min orbital pe- 
riod of 4U 1916—05 places it in the class of ultracompact 
LMXBs which re quire hydrogen-d eficient and/or degen- 
erate donors f^g.- lJoss et al.ll97^ . It has been proposed 
that the companion is a hydrogen-deficient but not yet 
degenerate star (Nelson et al. 1986). Observational ev- 
idence points to a helium-rich donor. The X-ray burst 
properties of 4U 1916—05 suggest that the accreting ma- 
terial is helium-rich and that the distance to 4U 1916—05 
is 8.9±1.3 kpc (Galloway et al. 2005, in prep.). A re- 
cent optical spectrum of 4U 1916—0 5 shows a feature at 
4540 A which may be due to He II IjNelemans fc .lonked 
|2QQ3). 

The non-dip spectrum can be fit with a power-law or 
power-law -I- blackbody model in the range 0.5-10 keV. 
Variations in the best-fit spectral parameters are corre- 
lated with the broad-band luminosity and t he position 
of 4U 1916—05 in its color-colo r diag ram (,Bloser et alJ 
|2Q03). In addition, lAsai et alJ l|200n|) detected a broad 
(FWHM = 0.7 keV) emission feature in the ASCA 
spectrum at 5.9 keV with an equivalent width of 87 
eV attributed to ¥e-K emission. Highly ionized iron 
absorption lines were found in the XMM spectrum of 
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4U 1916—05 with marginal detecti ons of other highly ion- 
ized features (e.g., Mg XII, S XVI: lBoirin et alJl2004ft . 

2. OBSERVATIONS AND DATA REDUCTION 

We observed 4U 1916-05 on 2004 August 07 for 50 
ks with Chandra using the High Energy Transmission 
Grating Spectrometer (HETGS) and the Advanced CCD 
Imaging Spectrometer (ACIS; iCanizares et all l2005[) . 
The HETGS carries two transmission gratings: the 
Medium Energy Gratings (MEGs) with a range of 2.5- 
31 A (0.4-5.0 keV) and the High Energy Gratings 
(HEGs) with a range of 1.2-15 A (0.8-10.0 keV). The 
HETGS spectra are imaged by ACIS, an array of six 
CCD detectors. The HETGS/ACIS combination pro- 
vides both an undispersed (zeroth order) image and dis- 
persed spectra from the gratings. The various orders 
overlap and are sorted using the intrinsic energy resolu- 
tion of the ACIS CCDs. The first-order MEG (HEG) 
spectrum has a spectral resolution of AA — 0.023 A 
(0.012 A). To reduce pileup, the observation used a sub- 
array of 512 rows, yielding a frametime of 1.7 s. The 
"level 1" event files were processed using the CIAO v3.2 
data analysis package^. The standard CIAO spectral re- 
duction procedure was performed. 

The combined MEG -f- HEG first order dispersed 
spectrum of 4U 1916—05 has an average count rate of 
10.6±1.3 cts s~^. During our observation, 4U 1916—05 
showed only mild dipping behavior as well as two X-ray 
bursts (see Figure ^). In the analysis presented here, 
we focus only on the persistent (non-dip) emission from 
4U 1916—05. We inspected the lightcurve at times of 
expect ed dipping, using the dip ephemeris of lChou et al.l 
(J200r). Of the 16 orbital cycles covered by our obser- 
vation, only six showed statistically significant dipping 
behavior. These dips were centered on a phase of w 0.15 
and had width of w 0.15 in phase. To produce our persis- 
tent (non-dip) spectrum, we excluded data from phases 
0.0-0.3 for all orbital cycles. In addition, we excluded 
data during the two X-ray bursts. We filtered the level 
2 event file, yielding a persistent (non-dip, non-burst) 
exposure time of 31 ks. We then extracted source and 
background spectra for both the plus and minus first or- 
der MEG and HEG. Detector response files (RMFs and 
ARFs) were created for the four spectra using the stan- 
dard CIAO tools. 

For bright sources , pileu p can be a problem for CCD 
detectors (see, e.g., iDavisi 12003). We checked the dis- 
persed spectra of 4U 1916—05 and found signs of pileup 
in the first order MEG spectrum. In dispersed spectra, 
pileup can affect only a limited wavelength range, par- 
ticularly where the effective area of the instrument is the 
highest. In our observation, pileup was present between 
2-10 A in the MEG plus and minus first orders. No 
pileup was found in the HEG spectra. 

The spe ctral analysis was performed using ISIS 
l)Houck fc D enicola 2000). The -1-1 and -1 order spec- 
tra were combined for the HEG and MEG respectively 
using the ISIS function combine_datasets. This func- 
tion is similar to summing the datasets, but more accu- 
rately takes into account the different responses for each 
dataset. We fit jointly the HEG spectra over the range 
1.6-11.5 A and MEG spectra over the range 10.0-14.0 A. 

^ http://asc.harvard.edu/ciao/ 



3. ANALYSIS AND RESULTS 

An initial inspection of the spectrum of 4U 1916—05 
revealed narrow absorption features at 12.1, 8.4, 6.2, 
and 4.7 A (see Figure |2), attributable to hydrogenic 
neon, magnesium, silicon, and sulfur. To determine the 
continuum model, we fit the spectrum ignoring regions 
around these lines as well as the Ve-K region at 1.8 A. 
We tested three different continuum models: power law, 
power law -f blackbody, and power law -I- disk black- 
body, all with interstellar absorption as described by 
the tbabs model using the abundances of I Wilms et alJ 
pOOO). Both the power law -t- blackbody and power law 
+ disk blackbody models give better fits to the data than 
the power law alone (see Table QJ. We take the power 
law -I- disk blackbody model as the continuum model, 
although we note that the choice of the power law -I- 
blackbody model would make little difference to our re- 
suhs. The absorbed 0.5-10 keV flux from 4U 1916-05 
was 8.0 X 10"^" erg cm~^ s~^. 

Our best-fit model is s i gnifica ntly different from 
that found b y iBoirin et aP Jlo^. But as noted by 
iBoirin et alJ (|2004D . the limited energy range covered 
by XMM and Chandra make it difficult to determine a 
unique solution to the continuum fit. One interesting dif- 
ference between our Chandra data a nd the XMM data, is 
the presence of an edge at 0.98 keV. IBoirin et al.l l|2004fl 
found that both the EPIC and RGS data were best-fit 
when an edge, with depth r = 0.11 ± 0.03, was included. 
In contrast, the addition of an edge at 0.98 keV does not 
improve the chi-squared value of the continuum fit for the 
Chandra data. We find an upper limit of r < 0.08 for the 
edge, barely consistent with the XMM result. The 0.5- 
10 keV luminosity of 4U 1916-05 was 9.0 x 10=^^ erg s-\ 
using a distance of 8.9 kpc, during the Chandra obser- 
vation. This is «2 times greater than found during the 
XMM observation l|Boirinet alJl2004r) . 

With the best-fit continuum model fixed, we fit the 
Hues from Ne X, Mg XII, Si XIV, S XVI, Fe XXV, and 
Fe XXVI with Gaussian models to determine the line 
positions, widths, and fluxes. In addition, we determined 
the upper limits to the flux from the helium-like ions of 
neon, magnesium, silicon, and sulfur. For the upper limit 
measurements, wc fixed the wavelength to those given 
in Bchar & Netzer (2002) and we fixed the line width 
(T — 0.005 A. We marginally detected the helium-like 
ion of sulfur. The neon lines were covered by the MEG 
spectrum, while the rest of the lines were located in the 
HEG spectral region. The best-fit parameters are given 
in Table 121 and the spectra and best-fit models are shown 
in Figures O and 01 

The line positions are consistent with the rest frame 
wavelengths of hydrogenic neon, magnesium, silicon, and 
sulfur, and hydrogenic and helium-like iron. This sug- 
gests that there is no bulk motion along our line of sight 
associated with the absorbing material. In addition, the 
measured line widths are comparable to, or less than, 
the instrument resolution (0.023 A for MEG, 0.012 A 
for HEG). Therefore, direct measurements of the velocity 
widths, or temperature, of the absorbing material are not 
possible with these data. The derived temperature limits 
range from < 4 x 10® K for Mg XII, to 1.2 ± 1.1 x lO^'' K 
for Fe XXVI. 

From the best fit equivalent widths (EW), we can esti- 
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mate the column density (A^z) of the ions using the rela- 
tionship between EW, Nz, and the transition oscillator 
strength fij 



A 



TTleC^ 



rA^zA/y = 8.85 x IQ-^^ jVzA/y , (1) 



where W\ is the EW in wavelength units and A, t he line 
wavelength, is given in cm units llSpitzgij ll978D . The 
equation above is true only when a line is unsaturated, 
on the linear part of the curve of growth. For saturated 
lines, the exact relationship between EW and Nz be- 
comes a more complicated function, and includes a de- 
pendence on the velocity of the ions. We used the os- 
cillator strengths given by Verner et al. (1996) for the 
hydrogenic lines, and those of Behar fc Netzer (2002) for 
the helium-like lines. Our results are given in Table[21 If 
the lines are saturated, the values represent a lower limit 
to the column densities for each ion. 

Com paring our results with those of iBoirin et al.l 
l)2004|) . we find that the iron line EWs are barely con- 
sistent within the errors. Our Fe XXVI EW is at the 
upper limit of the XMM result, while our Fe XXV EW is 
at the lower limit. It is possible that the greater luminos- 
ity of 4U 1916—05 during the Chandra observation has 
caused a real difference in the Fe XXV and Fe XXVI col- 
umn densities, compared to the lower luminosity XMM 
data, but the difference is not significant enough to con- 
firm this. The iron line widths found in the Chandra 
data are lower than the upper limits found in the XMM 
data owing to the higher resolution of the HETG com- 
pared to the EPIC pn. Our EW measurements for S XVI 
and Ne X are consistent with the values reported for the 
XMM data, but our Mg XII EW is significantly below 
the result of Boiri net al.l (|2004) . We note however that 
IBoirin et gJj (12004) concluded that their Mg XII detec- 
tion was only marginal, therefore the difference between 
the Chandra and XMM results may not indicate a real 
change between the two observations. As in the other X- 
ray dipper observations, we found no variation in the line 
parameters as a function of orbital phase (see Figure ISJ. 

4. DISCUSSION 

Our C/iarwira/HETGS observation of 4U 1916-05 re- 
vealed narrow, unresolved absorption lines in the persis- 
tent emission, attributable to hydrogenic neon, magne- 
sium, silicon, and sulfur, in addition to the previous iden- 
tified hydrogenic and helium-like iron absorption lines. 
This makes 4U 1916—05 only the second of the classical 
X-ray dipper systems to show narrow a bsorption lines 
from mid-Z elements (|Boirin et al.ll2004 and references 
therein). 

The properties of these lines is of particular interest 
for understanding the emission and absorption processes. 
With a 50 min orbital period, 4U 1916—05 is the shortest 
period X-ray dipper, and as such should have the small- 
est accretion disk. If the lines are associated with the 
accretion disk, they should have widths comparable to 
the expected Keplerian velocities in the disk. Assuming 
a 1. 4-2.0 Mp^ primary a nd a 0.1-0.15 Mq secondary (see, 
e.g. iNelson et al.lll986j) . we can estimate the outer disk 
velocity assuming that the disk fills 70% of the primary 
Roche lobe. We find that the outer disk velocity should 
be in the range of 1000-1300 km s^^, while the compan- 
ion velocity should be 740-840 km s^^. In all cases, the 



measured line widths are consistent with, or below, the 
instrument resolution. For Mg XII and Fe XXV, the up- 
per limits on the line widths are less than the estimated 
outer disk velocity. Additionally, we measure no shift in 
the wavelength of the lines, ruling out a strong outflow. 
The S XVI line shows a slight shift (-f 0.004 A) but given 
that this is the least significant of our line detections and 
the possibility of systematic errors, we do not feel it is a 
significant result. 

We suggest two possible explanations for these results. 
First, the absorber could be static, producing the narrow 
lines and lack of wavelengths shifts. This however seems 
unlikely given the large quantities of angular momentum 
present in the system. Instead, we propose that the line 
properties are a measure of the emission properties, par- 
ticularly the extent of the emission. In this case, we as- 
sume that the ionized absorber is associated with either 
the accretion disk or its atmosphere, as suggested by the 
implied cylindrical geometry (see e.g., Boirin et al. 2004). 
The rotation of the absorbing material is only measur- 
able when it has a component along the line of sight to 
the emitter. If the bulk velocity is perpendicular to the 
line of sight, no velocity effects will be found in the ab- 
sorption lines (excluding thermal or turbulent motions 
which we take as small compared to the Keplerian ve- 
locities in the disk). From the limits on the linewidths, 
we can calculate the maximum allowable radial extent of 
the emitter. We find that the X-ray emission region is 
< 32000 km, assuming that the absorber is located at 
the outer edge of the accretion disk. This represents the 
maximum upper limit on the extent. We note that if the 
absorber is located at a smaller disk radius, the upper 
limit will be smaller, roughly 1/4 of the absorber radius. 

Assuming that photoionization is the dominant ion- 
ization process for the absorbing material, we can esti- 
mate the ionization parameter, ^ = L/uer"^ (jTarter et alJ 
I1969|) . from the ratio of the helium-like to hydrogenic 
column densities for each atomic species. We assume 
that the helium-like and hydrogenic ions are present at 
the same, single-valued ionization parameter. This is a 
simple approximation to the expected true state of the 
system in which the ions exist in different relative pro- 
portions over a range of ionization parameters. 

We used XSTAR* to estimate the ion abundances as 
a function of ^ for an optically thin plasma, with con- 
stant number density and solar abundances. We note 
that while the companion in 4U 1916—05 is expected to 
be a hydrogen-deficient star, and solar abundances would 
therefore not apply, we found that varying the abun- 
dances made no significant difference in the helium-like 
to hydrogenic column density ratios for any element. The 
only XSTAR input that made a substantial difference in 
the inferred ionization parameters was the shape of the 
input spectrum. We tried three models, a power law with 
a = —1, a 10 keV bremsstrahlung, and a user defined 
model based on the continuum emission of 4U 1916—05. 

We compared the XSTAR results with the column den- 
sity ratios for neon, magnesium, silicon, sulfur, and iron 
as given in Tabled We note that these values are only 
valid if the lines are unsaturated. The limits for neon and 
magnesium are too high to provide any constraint on the 
ionization parameter. For silicon, we find log^ >1. 7-2.0 

* http://hcasarc.gsfc.nasa.gov/docs/softwarc/xstar/xstar.htmI 
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and for sulfur, log^ = 1.5-3.3, given the error in the 
measurement and the variation in ^ with spectral model 
(highest values for the power law model, lower for the 
bremsstrahlung and user defined models). The most con- 
straining result comes of course from the iron ratio, where 
we find log^ = 3.50±0.13 for the bremsstrahlung model, 
3.73±0.13 for the power law model, and 3.02±0.12 for 
the user defined model. This rough estimate would sug- 
gest that the iron absorption is found in a material with 
higher ionization than the lower Z elements. Our re- 
sult is consistent with modeling of the ionized absorber 
in MXB 1658—298 which required a distribution of ion- 



ization parameters l)Diaz Trigo et al.1 12006(1 . This is of 
course expected since in a material ionized enough to 
contain helium-like and hydrogenic iron, the lower Z el- 
ements would be completely stripped. 
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Table 1. Best-Fit Continuum Spectral Parameters'* 



Model 


A^H (lO^l cm-2 


r 


Alb 


kT (kcV) 


Norm'^ 


Xl/'^ 


PL 

PL+BB 

PL+DISKBB 


6.9±0.2 
5.6±0.7 
5.4±0.7 


1.519±0.019 

l.SliO.ll 

1.0±0.4 


12.2±0.3 

8.5±1.7 
5±3 


0.54±0.04 
0.91±0.11 


70±30 

15±5 


1.123/1187 
1.114/1185 
1.113/1185 



''All errors quoted at the 90%-confidence level. 

''Power-law normalization in units of 10~^ photons cm~^ s~^ keV~^. 
'^Thermal component normalization. For the blackbody model, norm 
(^in.km/^iOkpc)^ COS 6, where is the inclination angle of the disk. 



KJD^ 



lOkpc' 



for the disk blackbody model, norm 



Table 2. Best-Fit Line Parameters'' 



Wavelength (A) FWHM (A) FWHM (km s"!) EW (mA) Nz^ (10^^ cm'^) 



Ion 



NeX 


12.140±0.006 


0.035±0.019 


870±470 


-25±8 


4.6±1.5 


MgXII 


8.425±0.002 


<0.014 


<500 


-9±3 


3.4±1.1 


Si XIV 


6.185±0.002 


0.016±0.009 


800±460 


-11±3 


8±2 


SXVI 


4.736±0.003 


0.012±0.009 


750±600 


-7±3 


8±4 


FeXXV 


1.8519±0.0017 


<0.005 


<950 


-4.9±1.1 


20±5 


Fe XXVI 


1.7811±0.0015 


0.012±0.005 


1900±900 


-11±3 


90±30 


Upper Limits on Helium-like Ions 


NelX 


13.448 (fixed) 


0.012 (fixed) 


260 (fixed) 


>-9 


<0.9 


MgXI 


9.170 (fixed) 


0.012 (fixed) 


380 (fixed) 


> -6 


< 1.2 


Si XIII 


6.648 (fixed) 


0.012 (fixed) 


530 (fixed) 


> -3 


< 1.1 


SXIV 


5.039 (fixed) 


0.012 (fixed) 


700 (fixed) 


-4±3 


2.5±1.9 



"All errors quoted at the 90%-confidence level. 

''Assumes lines are unsaturated. Values represent a lower limit to the column density if lines are saturated. 
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Fig. 1. — Left panel, top: First order HEG-I-MEG 0.5-8.0 keV count rate over the Chandra observat ion of 4U 1916—05. The solid marks 
at the top of the figure indicate the expected dip times using the dip ephemeris offCJhou et alj 120011) . Left panel, bottom: The hardness 
ratio ([4—8 keV count rate]/[0.5— 1.5 keV count rate]) during the observation. As expected the hardness ratio increases during the X-ray 
dips (due to low energy absorption) and X-ray bursts. Right Panel: Normalized phase folded lightcurve of the first order HEG-I-MEG 
0.5—8.0 keV counts from 4U 1916—05. The times of X-ray bursts have been removed and the lightcurve has been corrected for exposure 
effects. The solid line from phases 0.3—1.0 indicates the data used to produce the persistent (non-dip) spectrum. 
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Fig. 2. — Flux spectra and best-fit model for the combined first order MEG+HEG Chandra spectrum of 4U 1916—05. The data and 
model arc binned to 0.010 A. The solid line is the full best-fit model including both continuum and line components. The dotted and dashed 
lines show the best-fit power law and disk blackbody continuum models, respectively. The absorption features from neon, magnesium, 
silicon, sulfur, and iron are labeled. 
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Fig. 3. — Flux spectra and best-fit models for the Ne X, Mg XII, Si XIV, and S XVI line regions. Note that the Ne X line is covered 
by the lower resolution MEG, while the other lines are covered by the HEG. The line-like feature at aill.97 A in the Ne X region is more 
narrow than the instrument resolution and is likely an instrumental artifact. 



Ionized Lines in 4U 1916—05 
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Fig. 4. — Flux spectrum and best-fit model for the iron line region. The Fe XXV line is found at 1.85 A and the Fe XXVI line is found 
at 1.78 A. 
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Fig. 5. — Comparison of the absorption line equivalent widths as a function of orbital phase. We plot the EWs for the full dataset 
(diamond), in addition to three phase groups: group A covering phases 0.30—0.53 (triangle), group B covering phases 0.53—0.76 (square), 
and group C covering phases 0.76-1.00 (x). No significant variation of the line properties are found. 



